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Abstract. Multiple parton interactions in relativistic heavy ion reactions result in transverse mo- 
mentum diffusion and medium induced non-Abelian energy loss of the hard probes traversing cold 
and hot nuclear matter. A systematic study of the interplay of nuclear effects on the pr >2 GeV 
inclusive hadron spectra demonstrates that the competition between nuclear shadowing, multiple 
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, scattering and jet quenching leads to distinctly different enhancement/suppression of moderate and 

high-pr hadron production in d+Au and Am + Am collisions at RHIC. The associated increase of 
di-jet acoplanarity, measured via the broadening of the back-to-back di-hadron correlation function, 
provides an additional experimental tool to test the difference in the dynamical properties of the 
^ ' media created in such reactions. 
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. Particle production from a single hard scattering with momentum exchange much 

larger than 1/fm is localized in space-time. It is multiple parton scattering before or 
after the hard collision that is sensitive to the properties of the nuclear matter [1, 2]. By 
comparing the high-pj- observables in p + p, p +A and A+A reactions, we are able 
^ ! to study the strong interaction dynamics of QCD in the vacuum, cold nuclear matter 

■ and hot dense medium of quarks and gluons, respectively. So far the first two integral 

moments / dz z'^p (z) of the matter density in the interaction region can be deduced from 
^ '. experimental measurements since they are related to the broadening [3, 4] and energy 

H ' loss [5, 6] of a fast parton traversing nuclear matter. From [3, 5] around midrapidity: 
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In Eq.(l) the factor 2 comes from 2D diffusion, ^ ~ ^(1) and p^ is the effective 
gluon density. For the 1-i-lD Bjorken expansion scenario A^ is the transverse area of 
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FIGURE 1. Left panel: predicted suppression for 71" and /i+ +/z [1] in Am + Am compared to PHENIX 
and STAR data [7]. Similar quenching is found by BRAHMS and PHOBOS [7]. Right panel: predicted 
small enhancement of neutral pions in d +Au [8], data is from PHENIX. Left bottom panel: a test of a 
suggested interpretation of high /^^-hadron suppression as result of initial state wavefunction modification 
(RdAu ~ VRZIm,)- Data is from BRAHMS, PHOBOS and STAR [8]. 



the interaction region, Tq is the initial equilibration time and dN^/dy is the effective 
gluon rapidity density. In Eq.(2) the dominant logarithmically enhanced contribution to 
mean energy loss computed in the GLV approach [5] is shown. For further discussion 
on medium induced radiative energy loss in QCD see [5, 6]. 

Initial state parton broadening, nuclear shadowing and jet energy loss are incorporated 
in the lowest order pQCD hadron production formalism as in [1]. The interplay of 
dynamical nuclear effects can be studied through the nuclear modification ratio 
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Rab{pt) = -, — 7i — / T-p^ , about impact parameter b in A +5 . (3) 

dyd'^^T dyd'^^T 

Figure 1 compares the predicted [1] approximately constant suppression of n'^ and 
h+ + h- in V5 = 200 AGeV Au+Au collisions at RHIC to PHENIX and STAR 

data [7]. The overall quenching magnitude and its centrality dependence are set by 
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( (L) /A^)dN^ /dyoc Npan, dN^ /dy=\\50. The shape of RauAu is a result of the interplay 
of all three nuclear effects. The full numerical calculation takes into account the dynam- 
ical Bjorken expansion of the medium, finite kinematic bounds, higher order opacity 
corrections and approximates multiple gluon emission by a Poisson distribution [1, 5]. 
In the right panel the Cronin enhancement, resulting from initial state parton broadening, 
vac ~ il^^ / ^g)eff{L) [1, 3] is seen to compare qualitatively to the shape 
of the PHENIX measurement [8] in d +Au. The calculations in Fig. 1 use the value 
{ll^/Xg)eff = 2 X 0.14 GeV^/fm for the cold nuclear matter transport coefficient con- 
strained from existing Cronin data [1], although somewhat smaller scattering strength 
may be favored by PHENIX data. Larger enhancement of /i+ -|- production, consis- 



tent with results form low energy p+A measurements, is also shown [8]. The lower 
right panel rules out the scenario for the initial wavefunction origin of moderate and 
high-pT hadron suppression (see left panel of Fig. 1) since in this case RdAu ~ V^AuAu- 
For further discussion on the Cronin effect see [1, 9]. 



Broadening of the away-side di-hadron correlation function 

The total vacuum+nuclear induced broadening for the two partons in a plane perpen- 
dicular to the collision axis in p+A (A+ A) reads [3]: 

(4) = (4) vac + .2^"\ (^) {L)lS + 2jets(^^) (^) {L)fS ■ (4) 

[^jets) \^ Jeff V2/ projection \ ^ / eff 

A typical range for the cold nuclear matter transport coefficient for gluons is given by 
{H^/Xg)eff,is^FS = 2 X 0.1 GeV^/fm - 2 x 0.15 GeV^/fm. For final state scattering in 
a 1+lD Bjorken expanding quark-gluon plasma the broadening {{ix'^ / X)eff{L)ps) can 
be evaluated from Eq.(l) with dN^/dy = 1150, consistent with the inclusive hadron 
suppression pattern in Fig 1. 
Figure 2 shows two measures of the predicted increase in di-jet acoplanarity for 

minimum bias d +Au and central Au+Au reactions [3]: (|kry|) = ^J {k^)iparton/^, 
I parton — (kr)/2) and the away- side width Opar of the di-hadron correlation func- 
tion C(A(^) = N^^'''^{A^)/Nl'^t^\ C{A^) is approximated by near-side and far-side 
Gaussians for a symmetric pj^ fsi pf- case and the vacuum widths are taken from 
PHENIX [10]. In the right panel of Fig. 2 di-hadron correlations in J + Am are shown to 
be qualitatively similar to the p+p case and in agreement with STAR measurements [8]. 
In Au+Au reactions at RHIC di-jet acoplanarity is noticeably larger, but this effect alone 
does not lead to the reported disappearance of the back-to-back correlations [10]. To 
first approximation the coefficient of the away-side Gaussian (the area under C{A^), 
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A(j) > n/2), is determined by jet energy loss and given by Raa °^ Np^rt- Broadening with 
and without away-side quenching is shown the bottom right panel of Fig. 2. Combined 
<i-|-Au and Am -I- Am experimental data in Fig. 2 also rule out the existence of monojets 
at RHIC. For further discussion on di-hadron correlations see [3, 11]. 

In summary, evidence from jet tomography [1,2], relativistic hydrodynamics [12] and 
parton cascade models [13] is in strong support of the creation of a deconfined phase 
of QCD at RHIC with initial energy density ~ 20 GeV/fm^, more than 100 times the 
1 /7 GeV/fm^ density of cold nuclear matter. 
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FIGURE 2. Left panel: predicted enhancement of (|k7-y|) and a far in minimum bias d + Au and 
central Au + Am reactions at RHIC from pj- -diffusion [3]. Preliminary p + p data is from PHENIX [10]. 
Right panel: the broadening of the far-side di-hadron correlation function in central d +Au and Au +Au 
compared to scaled (xlO) STAR data [8]. In the bottom right panel the broadening with and without 
suppression, approximately given by R^a from Fig. 1, are shown. 
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